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Abstract 

We analyze the Eckhaus instability of plane waves in the one-dimensional com- 
plex Ginzburg-Landau equation (CGLE) and describe the nonlinear effects arising 
in the Eckhaus unstable regime. Modulated amplitude waves (MAWs) are quasi- 
periodic solutions of the CGLE that emerge near the Eckhaus instability of plane 
waves and cease to exist due to saddle-node bifurcations (SN). These MAWs can be 
characterized by their average phase gradient v and by the spatial period P of the 
periodic amplitude modulation. A numerical bifurcation analysis reveals the exis- 
tence and stability properties of MAWs with arbitrary v and P. MAWs are found to 
be stable for large enough v and intermediate values of P. For different parameter 
values they are unstable to splitting and attractive interaction between subsequent 
extrema of the amplitude. Defects form from perturbed plane waves for parameter 
values above the SN of the corresponding MAWs. The break-down of phase chaos 
with average phase gradient v ^ ^ ("wound-up phase chaos") is thus related to 
these SNs. A lower bound for the break-down of wound-up phase chaos is given by 
the necessary presence of SNs and an upper bound by the absence of the splitting 
instability of MAWs. 
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Fig. 1. Phase diagram of the one-dimensional CGLE. The dotted curve indicates the 
Benjamin-Feir-Newell hne. Plane waves undergo an Eckhaus instability at values of 
ci, C3 below this curve depending on their wavenumber. Above the dashed curve the 
Eckhaus instability is supercritical whereas it is subcritical below [4]. MAWs and 
wound-up phase chaos with 1/ > can be observed between the dashed and the full 
curve. Defect chaos can occur only above the full curve [27,28] which denotes the 
saddle-node bifurcation of MAWs with v = and P 00. The vertical dot-dashed 
line indicates the cut of the parameter space at ci = 3.5 studied in this paper. 

1 Introduction 



The emergence of chaotic behaviour from ordered states in spatially extended 
systems has been the subject of many recent experimental and theoretical in- 
vestigations [1,2]. Nonetheless, the mechanisms leading from stationary regimes 
to chaotic (or spatially irregular) phases still pose many challenging questions. 
One of the most studied instabiUties in extended oscillatory systems is the 
Eckhaus instabihty of plane waves [3] . 

The occurrence of this instability has been experimentally observed in many 
quasi one-dimensional systems like the oscillatory instability of a Rayleigh- 
Benard convection pattern [4], hydrothermal waves [5-8], heated wire convec- 
tion [9], sidewall convection [10], the Taylor-Dean system [11] and internal 
waves excited by the Marangoni effect [12]. The Eckhaus instability also plays 
an important role in the radial dynamics of spiral waves in the Belousov- 
Zhabotinsky reaction [13]. 

The complex Ginzburg-Landau equation (CGLE) [2,14] is the appropriate am- 
plitude equation to describe the slow dynamics near a supercritical transition 
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to unidirectional traveling waves. In one spatial dimension, the CGLE reads: 

dtA = A+{l + tci)dlA - (1 - ic3)\A\^A , (1) 

where Ci and C3 are real coefficients and the field A = A{x, t) = \A{x, t)\e^'^^^'^^ 
has complex values. As exact solutions the CGLE admits plane waves of the 
form Aq{x,t) — age*^*^~'^«*\ where q indicates the wavenumber, = — q^ 
and Uq = -C3 + g^(ci + C3). 

A linear stability analysis [15] of these solutions can be performed by consider- 
ing the perturbed solution Aq{x, t) — (a^ + 5a)e*(^^"''^«*\ where 5a oc q^^^q'^^^)*. 
The growth rates associated to the complex perturbation 5a is 

±^J{l + ci)(l - g2)2 _ - 2iqk - c^{l - q^)f . (2) 

The plane waves become linearly unstable to long wavelength perturbations 
{k 0) for q = qE = sj {I - Cic^)/{2{1 + ci) + 1 - C1C3). This limit is called 
Eckhaus instability [3]. Above the Benjamin-Feir- Newell (BFN) line 1 — C1C3 = 
0, all plane waves are unstable to homogeneous perturbations. For a given q > 
qE, the corresponding plane wave is linearly unstable against perturbations 
with wavenumbers k inside the interval < \k\ < kc- kc increases for increasing 
values of the parameters ci and C3 [16]. 

As noticed in [17,18], the Eckhaus instability is a convcctivc instability. Thus, 
it is relevant in the systems with periodic boundary conditions considered 
here, while it would be suppressed in fixed boundary conditions , e. g. zero- 
flux or Dirichlet boundary conditions. In the latter geometries, the absolute 
instability of the plane waves has to be computed. It occurs for sufficiently 
large q or/and ci,C3 values inside the Eckhaus unstable range. 

Another interesting aspect of the Eckhaus instability is found when the non- 
linearities of the CGLE are taken into account. A weakly nonlinear analysis 
[4] revealed that for 

cl{l - 6cl) + ci{2cl + I6C3) - (8 + c^) > (3) 

the Eckhaus instability becomes supercritical, i.e. the instabilities are satu- 
rated and the emerging quasi-periodic solutions {resp. modulated amplitude 
waves) coexist with the unstable plane waves. Numerical simulations [4,19-22] 
provided examples of such stable modulated amplitude waves (MAWs) . Stable 
MAWs have also been observed in experiments on surface-tension-driven hy- 
drothermal waves [6] as well as on the Taylor-Dean system [11] and on internal 
waves excited by the Marangoni effect [12]. 
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Fig. 2. Maximum (dots) conserved average phase gradient i^mCcs) for fixed Ci = 3.5 
obtained from numerical simulations (with system size L = 1,024 — 2,048 and 
integration times t ^ 10^) for 50-70 different initial conditions (noise added to 
plane wave with wavenumber q = v). For v < vm no defects were present while 
above vm at least one initial condition caused defects. The full curve denotes the 
Eckhaus instability of plane waves that converges to the BFN line at C3 = 1/ci for 
1/ = 0. For V above the diamonds regular states were observed after a transient phase 
chaotic dynamics but below the diamonds most initial conditions led to persistent 
spatio-temporal chaos [21]. Li denotes the lower bound for the occurrence of defect 
chaos in the thermodynamic limit as calculated in [27,28]. 

In addition, the CGLE exhibits two qualitatively different spatiotemporal 
chaotic states known as phase chaos (when the modulus of the field \A\ is 
bounded away from zero) and defect chaos (when the phase of A displays 
singularities where \A\ =0) [23-26]. The subclass of MAWs with zero average 
phase gradient is important for understanding the transition from phase to 
defect chaos (see full curve in Fig. 1) [27,28]. In the phase chaos regime states 
with nonzero average phase gradient v have a dynamics quite different from 
that at z/ ~ 0. In particular, these states can be either chaotic or regular de- 
pending on the initial conditions and on the parameters Ci,C3 and v. In this 
paper we will focus on MAWs with v ^ Q and on the dynamical regime asso- 
ciated to them, that is referred to as "wound-up" phase chaos [20]. It will be 
shown that MAWs and wound-up phase chaos exist between the dashed and 
the full curve in Fig. 1. 

In Section 2 the analysis of MAWs introduced in [27,28] is extended to arbi- 
trary values of the average phase gradients of the field. The two parameter 
family of MAW solutions is parametrized by the spatial period P of the mod- 
ulation and by the average phase gradient 
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For plane wave solutions, v equals the wavenumber q. In analogy, the phase 
gradient is often called "local wavenumber" . A linear stability analysis will 
show that MAWs with i/ 7^ can be stable even in infinitely large systems. 
In contrast, MAWs with v — ^ are always unstable in systems of large length 
L ^ P. For systems with periodic boundary conditions the average phase 
gradient of the whole system can only be changed, if a space-time defect occurs 
: drops to zero and (fx locally diverges at a defect. Persistent phase 

chaos with conserved u < vm 7^ has been observed in numerical simulations 
of the CGLE (1) [20,21]. The maximum conserved average phase gradient 
vm decreases as function of the coefficients Ci,C3 [20,21] and vanishes at the 
apparent transition from phase to defect chaos, vm was therefore suggested 
[21] as an order parameter for this transition. We extended the numerical 
determination of vm towards smaller C3 and report the corresponding data in 
Fig. 2. 

In Section 3, a nonlinear analysis of the Eckhaus instability allows estimates 
for which parameter values defects occur. Lower and upper bounds for the 
limit Vm of wound-up phase chaos are derived from the existence and stability 
properties of the MAWs. For increasing v the degree of chaoticity associated 
with the wound-up phase chaos decreases [21]. For large enough u the dynam- 
ics can even become regular and stable quasi-periodic MAWs appear. The 
diamonds in Fig. 2 indicate this stability limit for numerical simulations with 
fixed ci = 3.5. The analysis in Subsection 2.3 will clarify this observation. 

The large number of parameters (ci, C3, z/, P) calls for restrictions. We limit our 
analysis to fixed ci = 3.5 since most previous numerical work has been done at 
this value [21 ,26] . The results will be presented as projections of the P direction 
onto the (cs, u) plane as well as in cuts through the parameter space spanned 
by C3,z/ and P. Additional investigations of the existence domains of MAWs 
revealed qualitatively similar results for fixed Ci = 0.4, 1.2, 2.1 and 5 and vari- 
able C3 as well as for fixed C3 = 0.83 and variable Ci. Two of these choices were 
studied by numerical simulations in [20]. A similarity transformation maps 
coherent structures onto each other along curves (ci -|- C3)/(l — C1C3) = const 
in coefficient space [14]. The parameters v, P, lo, v are transformed accordingly. 
One can thereby extend the results presented here to other values of the coef- 
ficient Ci. Section 4 discusses possible observations of MAWs in experimental 
systems. Finally, Section 5 summarizes the main results. 



2 Existence and stability of modulated amplitude waves 
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2.1 Coherent structure approach 

In analogy with the hnear analysis of the Eckhaus instability of plane waves 
we make the following ansatz for saturated modulations 

>l(x,i) = a(^)e^^We^(^^-'^^) (5) 

and rewrite it as 

yl(a;,^) = a(z)e^'^(^)e^^^ (6) 

where a and are real-valued functions of z := x — vt and (f){z) = (j){z) + 
qz,uj = qv — uj. Here a{z) and (f){z) represent coherent structures [29]. Coherent 
structures have been studied extensively [21,27-30] and play an important role 
in various regimes of the CGLE [4,19-21,27-31]. 

Substitution of ansatz (6) into the CGLE (1) yields the set of three coupled 
nonlinear ordinary differential equations (ODEs) 

b^^i/j'^a - -f-'^[{l + ciu})a + v{b + ciipa) - (1 - CiC3)a^] (7) 
■0^ = -26'0/a + 7~'^[ci —u; + v{cib/a — ip) — (ci + cs)a^] 

where b: = az,ip:=4>z and ■j:=l + c1 [32]. The continuation software AUT097 
[33] is used to compute the periodic orbits of the ODEs (7) that correspond 
to spatially periodic functions a{z), (f){z). In order to choose a unique solution 
from the continuous two-parameter family of periodic orbits we set the system 
size L equal to the period P of the periodic orbit and fix its average phase 
gradient hy v = Jq ijjdz . 

The continuation procedure starts from a fixed point (a, b, ■0) = — q'^, 0, q) 
that corresponds to a plane wave solution. Varying C3, a Hopf (HB) bifurca- 
tion (filled square in Fig. 3) is detected in the ODEs where the mode with 
the smallest possible wave number knB = ^n/P destabilizes the plane wave. 
Continuing the resulting branch of MAWs the free parameters uj and v are 
adjusted by the continuation algorithm. The continuation follows a unique 
branch of MAWs with u — q and P — L. Fig. 3 shows examples of resulting 
bifurcation diagrams. 
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Fig. 3. (a) Example of a bifurcation diagram showing the maximum of the modulus 
for MAWs with v = 0.25, ci = 3.5, P = li^jv = 25.13. The plane wave that is stable 
(unstable) against modes of wavelength P is represented by the thin full (dashed) 
line. The stable lower branch (unstable upper branch) of MAWs is denoted by the 
thick full (dashed) curve. HB denotes the Hopf bifurcation (square) of the plane 
wave solution whereas SN stands for the saddle-node bifurcation (triangle) that 
limits the existence of MAWs. Spatial portraits of (b) the modulus and (c) the 
phase gradient are shown for a choice of solutions. The dotted line represents the 
plane wave whereas thin full and thick full curves give MAWs at locations labelled 
by A and SN in (a). The dashed curve denotes the saddle- type upper branch solution 
at C in (a), (d) shows the oscillation frequency lo = qv — uj and (e) the velocity v 
versus C3. In (c) the dotted line denotes the group velocity (k = 0) and the line 
below gives the velocity Vc corresponding to the mode with finite wavelength P 
[16]. 



2.2 Existence limits of MAWs 



Upon increasing of C3 amplitude modulations grow and develop a localized 
depression where (px has a maximum (see Fig. 3b, c). As for u — 0, these 

MAWs are called the lower branch in contrast to the coexisting upper branch 
MAWs. The upper branch MAWs arc always unstable, while the MAWs of 
the lower branch can be stable in appropriate parameter regions. Examples of 
these lower branch MAWs have been obtained by numerical simulations earlier 
[4,19-21]. They have been analyzed in detail in [22]. Numerical simulations 
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Fig. 4. Bifurcation diagram with exclusively upper branch MAWs for 
V = 0.68, i-* = 27r/z/ « 9.24. The solid line indicates stable and the dashed curves 
unstable solutions. The Hopf bifurcations are both subcritical. 

can neither uncover unstable upper branch MAWs nor elucidate the limits 
of existence of MAWs. The bifurcation analysis presented here reveals that 
upper and lower branch meet and terminate in a saddle-node (SN) bifurcation 
(filled triangle in Fig. 3 a.d.c ). Due to the SN bifurcation the upper branch 
MAWs always have at least one unstable eigenmode, see also [30,31]. The 
upper branch continues to negative C3 and there connects to another instability 
of the plane wave with identical v and P. In the following we will concentrate 
on the lower branch MAWs. 

For large v and small P, the Hopf bifurcation is no longer supercritical and an 
unstable branch emerges directly from the plane wave. This is in agreement 
with analytical predictions [4]. Fig. 4 shows an example which also includes 
the second HB at negative C3. For v — ^ the MAWs emerge stationary [27,28] 
and acquire t; 7^ above a subsequent drift pitchfork bifurcation [34]. In the 
present case z/ 7^ the plane wave already breaks the reflection symmetry, the 
initial MAW has a nonzero velocity and the drift pitchfork (DP) bifurcation 
(filled diamond) is unfolded. See Fig. 5 for an example at fixed C3 = 2. The 
branch emerging at the HB in Fig. 5b represents the MAWs as discussed 
above. The second branch in Fig. 5b emerges at the period doubling (PD) 
bifurcation (open square) of MAWs with half the period. It always has unstable 
eigenmodes that drive the dynamics away from it to the coexisting MAWs of 
shorter period. Therefore this (upper) branch plays no essential role and is 
not treated further. 



2.2.1 Infinite system size 

We have analyzed the existence of lower branch MAWs in the entire parameter 
space (c3, u, P) at fixed Ci = 3.5. The system size is assumed infinitely large 
in order to allow for arbitrary periods P of MAWs. Fig. 6 shows examples of 
existence domains for P = 15, P = 30 and P — > 00. We find that both HB 
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Fig. 5. Bifurcation diagrams showing the velocity v versus ci. (a) Branches with 
V emerge at the drift pitchfork bifurcation (DP) for C3 = 2, P = 25, = 0. (b) 
The bifurcation is unfolded for v ^ 0, here C3 = 2, P = 25, 1/ = 0.01. An equivalent 
pair of branches exists ioi v ^ —u and v ^ —v. 

and SN shift to larger C3 as the period P is decreased. The same behavior has 
already been observed in the special case z/ = [27,28]. The dotted curve in 
Fig. 6 indicates the "envelope" of all SN bifurcations for MAWs of any period 
and therefore is the upper boundary for the existence domain of the MAWs. 



2.2.2 Medium system size 

Experimental setups and numerical simulations are restricted to finite system 
size L. Often periodic boundary conditions (corresponding to an annular ge- 
ometry) are used in order to study bulk effects of extended systems and to 
minimize boundary effects. The periodic boundary conditions also restrict pos- 
sible modes of perturbations. As described by Eq. (2) the instability threshold 
of plane wave solutions depends on the wavenumber k of the perturbation. 
Since in the studied range of coefficients the Eckhaus instabihty is a long- 
wavelength instability the plane waves will be stabilized in small systems. 
The instability threshold is shifted to larger values of the coefficients Ci,C3 
and can be computed from Eq. (2) setting k = 27t/L. 

Clearly the selection of perturbations by periodic boundary conditions also 
restricts possible MAWs. Their average phase gradient u and the period P 
have to be consistent wdth the system size and this renders the two-parameter 
family of MAWs discrete. It is thus convenient to parametrize MAWs by the 
average phase gradient u and the ratio n of wavelength 

P 

n := —- . (8) 

The ratio n takes values of integer fractions where the nominator counts the 
number of underlying wave length and the denominator the number of humps 
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Fig. 6. Existence domains of MAWs with period P projected onto the [03, u) pa- 
rameter plane. Thin curves denote the Eckhaus instabihty and HB which occur 
supercritical (full curve) or subcritical (dashed) depending on ly and P. The thick 
curves give the SN for selected periods P. Three examples of existence domains for 
P = 15 (right, dark shaded domain), P = 30 (middle, empty), P ^ 00 (left, light 
shaded) are shown. The superposition of all existence domains is bounded by the 
dotted curve. 

of the modulation. Hence this quantity is easily accessible in experiments. 
The existence domains of MAWs with respective n are presented in the (03, i^) 
parameter plane in Fig. 7. 

2.2.3 Small system size 

Here we focus on the extreme case. The shortest possible system with pe- 
riodic boundary conditions only contains one wavelength of the plane wave, 
consequently its length L is given hy L — 2t:/i/. In [21] the quantity vu was de- 
termined in analogy to vm for large systems, vu denotes the largest v for which 
none of the random initial conditions (different realizations of noise added to a 
plane wave) produced a defect. In the following these data (symbols in Fig. 8) 
are compared to the existence domains of MAWs. 

Within the light shaded area in Fig. 8 plane wave solutions with wavenumbcr 
V are stable in the short system. The stability area extends over the phase 
chaos and into the defect chaos region. We stress this result because from 
the experimental observation of stable plane waves one cannot necessarily 
infer that the dynamics of the system may be reproduced by the CGLE with 
coefficients Ci, C3 in the Benjamin- Feir stable region. The dashed curve denotes 
a subcritical instability. Only unstable upper branch MAWs exist to the left 
of this curve. For smaller u the instability again turns supercritical and stable 
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Fig. 7. Existence domains of lower brancli MAWs are denoted by shaded areas for 
(a) n = 20, (b) n = 4, (c) n = 2, (d) n = 1, (e) n = 1/2 and (f) n = 1/4. They are 
Umited by HB (sohd curve) at small C3 and by the SN (dashed curve) at large C3. 
Spatial profiles of coexisting MAWs at v = 0.05, C3 = 0.5 are shown for (g) n = 1, 
(h) n = 1/2 and (i) n = 1/4, corresponding to dots in (d-f). 



lower branch MAWs exist inside the dark shaded region. The thick full curve 
gives the SN bifurcation for MAWs with P = L. The thin curves show the 
respective limits of MAWs with shorter period. Defects are expected beyond 
the SN [27,28] and the subcritical instabihty which well reproduces the data 
from numerical simulations [21] except at small i^. Simulations with v < vjj 
resulted in modulations with a single hump (squares) or with two (triangles) 
or three (stars) humps of different size [21]. The latter two are observed above 
the SN of MAWs with P = L. Here the initial conditions select MAWs with 
shorter period which only coexist at small v. The SN with P — L nevertheless 
gives a lower bound for the formation of defects. 



2.3 Instabilities of MAWs 



In contrast to the case of MAWs with zero average phase gradients some 
MAWs with non-zero i> are stable even in very large systems. A linear stabil- 
ity analysis of MAWs as in [28] yields the spectrum of eigenvalues as shown in 
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Fig. 8. For short system size L = 27r/z/ plane waves are stable for parameter choices 
inside the light shaded area. MAWs with a single hump (P = L) exist inside the 
dark shaded area bounded by the supercritical HB to the left and the SN to the 
right. Thin curves give the limits of MAWs with two humps P = L/2 (dashed) 
and three humps P = L/3 (dotted). See the legend for the different cases. Sym- 
bols denote maximal u = uu that did not create defects but resulted in stable 
asymptotic states in simulations of the short system. Plane waves and single MAWs 
(squares), multi-hump MAWs with 2 humps (triangles) and 3 humps (stars) have 
been observed at v = vjj. Data ujj are taken from [21]. 

Fig. 9 for a typical example. From Fig. 9 we conclude that for this example the 
entire spectrum in the infinite system will be confined to the left half-plane. 
Thus MAWs should be found in experiments, that can be well described by 
the CGLE for appropriate control parameters. In this section we present a de- 
tailed study of the stability properties of MAWs. MAWs with a single hump 
per period P will be called "single MAWs". Their existence domains were 
studied in the previous section. However, the effective interaction between 
adjacent periods of a single MAW can be repulsive or attractive (see Figs. 
11 and 12). Period doubhng (PD) bifurcations (open squares) occur at the 
transitions from repulsive to attractive interaction [36]. There, new branches 
of MAWs with longer period but many humps per period emerge from the 
primary branch of single MAWs. We will call these solutions "multi-hump 
MAWs". In their profile some humps gain more space and others are com- 
pressed in an alternating fashion. The new branches extend to larger ca than 
the corresponding single MAWs. Fig. 10 shows how these branches arrange 
in a system with 4 interacting humps (L = 4 * P). As long as the PD bifur- 
cations are supercritical, the multi-hump MAWs are stable. They represent 
the saturated solution for attractive interaction between subsequent modula- 
tions. For large systems a whole sequence of period doubling bifurcations will 
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Fig. 9. (a) Spectrum of eigenvalues A of a stable lower branch MAW. Parameters 
are C3 = 0.4, = 0.184, P = 2'k/u. (b) Blow-up of the leading part of the spectrum. 
The dots correspond to system size L = 100 * P = 3415 and have been calculated 
using the Bloch method [28,35]. 

lead to multi-hump MAWs with an overall period equal to the system size. 

Hence they appear as an erratic spatial sequence of humps and depressions. 
This spatial sequence propagates in a coherent fashion. We named these pat- 
terns multi-hump MAWs to emphasize the connection among the coherent 
structures. 

Examples of these stable aperiodic patterns were already observed in numerical 
simulations. R. Montague et al. [20] denote this behavior as "frozen phase 
turbulence" while A. Torcini et al. [21] use the term "solutions of type 

The observed coexistence of a large number of stable multi-hump MAWs re- 
sults in a strong dependence of the final state on the initial conditions of the 
numerical simulation. Although each regular final configuration must be con- 
sistent with a particular single or multi-hump MAW it is difficult to predict 
how the selected final patterns depend on the initial conditions. 

MAWs with large period P undergo a "splitting" instability as in the limit case 
= [28] . Roughly, the spatial profiles of these MAWs consist of a localized 
hump and a plane wave part. Since the extended plane wave is linearly unstable 
the splitting instability is reminiscent of the Eckhaus instability. It creates 
more humps on the plateau of the unstable MAW and reduces the period P 
of MAWs on average. 

Figs. 11 and 12 represent cuts through the parameter space at fixed v = 0.25 
and C3 = 0.5, respectively. They show the typical arrangement of stable and 
unstable parameter regions for single MAWs. Other examined cuts for C3 — 
0.1, 0.2, 0.3, 0.4, 0.6 and 0.7 qualitatively show the same order. 

The cut through parameter space C3, u, P at u = 0.25 is shown in Fig. 11. The 
HB (dashed curve in the figure) approaches the Eckhaus instability for P 00 
as the lower bound of the existence domain. From above the domain is limited 
by the SN (sohd curve). For small P (large C3) the HB is subcritical and only 



13 



20 

max |A|^ 
15 

10 



0.5 



(a) 


|(d) 




/^multi-hump h 


single \ 


f] MAWs 1 


MAWs I 


1 // 


(b) y l/j 

/ / / 






plane wave 


1 . 1 . 1 . 1 . 



0.6 



0.7 0.8 

c. 



0.9 



0.5 



|A| 



(b) 

A ,A 


V \ 


(c) 




(d) 

-A A A (V^-^. 

I.I.I 



20 40 60 
z 



Fig. 10. (a) Bifurcation diagram for ly = 0.184, P = ir/iy and L = 4 * P = 68.3 
hence 4 pulse-like modulations interact. The maximum of the amplitude gradient is 
plotted since the interaction causes pulse shifts and the amplitude of single humps 
changes little. Thicker lines correspond to smaller overall period of the modulation. 
Typical solutions are shown in (b)-(d) as indicated by arrows in (a). 

unstable upper branch MAWs exist to the left. In the infinite system MAWs 
are found to be linearly stable for a broad range of parameters (dark shaded 
area). At low P the interaction instability occurs (white area) whereas at large 
P the long plateau of the MAW is unstable to splitting (hght shaded area). 
For C3 < 0.45 most random initial conditions will evolve to stable MAWs. 



A cut perpendicular to the previous one is shown in Fig. 12. Curves and 
shadings have the same meaning as discussed above. Starting from random 
initial conditions at > 0.1 a transient may again lead to a stable MAW 
with local periods P inside the stable windows. At lower u < 0.1 the prob- 
ability of approaching a stable configuration decreases since only a third of 
the previous stable P intervals remains. Below u = 0.02 no stable state can 
be prepared at all. Instead one observes wound-up phase chaos with an as- 
sociated maximal Lyapunov exponent that increases for decreasing v. As in 
the limit case 1/ — (phase chaos), the dynamics is driven by the attractive 
interaction and annihilation of localized modulations in competition with the 
splitting instability that produces new peaks in the modulations. In partic- 
ular for decreasing u, the splitting instability extends to shorter periods P 
and significantly overlaps with the interaction instability. With the above ar- 
guments many results obtained by numerical simulations of the CGLE can 
be well interpreted. In particular it has been observed in [21] for the same 
choice of parameters (ci = 3.5 and C3 = 0.5) that the maximal Lyapunov ex- 
ponent (averaged over many different initial conditions) is positive for 1/ = 
and decreases monotonously towards zero for increasing i/. Above u — 0.09 no 



14 




Fig. 11. Stability domain (dark area) of single MAWs for v = 0.25, L — oo. MAWs 
exist between supercritical HB (dashed curve) and SN (full curve). The tick marks 
at the right frame give the asymptotic values for P ^ oo. The dot-dashed curve 
denotes the subcritical HB. MAWs are unstable to splitting within the light shaded 
domain at large P. Within the white domain at small P single MAWs are unstable 
to interaction. 



chaotic solutions have been observed. 



3 Defect formation in wound-up phcise chaos 



In this section the formation of defects and the resulting change of the average 
phase gradient are studied. For v 0, the scenario of defect formation past the 
saddle-node bifurcation of the relevant MAW is analogous to the previously 
studied case u = [27,28], see Subsection 3.1. In particular, the dependence 
of the final selected average phase gradient Uf on the initial value z/j found in 
numerical simulations [20] can be interpreted. Subsection 3.2 is then devoted 
to the limit um of wound-up phase chaos. For a certain range of parameters the 
limit I'M is reproduced by means of the stability properties of MAWs at the 
saddlc-node bifurcation. These arguments work well for defect creation with 
u > 0.1. At smaller values of u various instabilities (splitting and interaction) 
of MAWs compete and a general statement is more difficult, compare also [28]. 
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Fig. 12. Stability domains (dark areas) of single MAWs for C3 = 0.5, L 00. Curves 
and instability domains (splitting=light shaded, interaction=white) have the same 
meaning as in Fig. 11. 

3.1 Beyond the saddle-node bifurcation 



The role of the SN bifurcation for the dynamics has been studied in [27,28] for 
the limit case 1/ = 0. For 1/ ^ Owe find similar behavior. Fig. 13 gives examples 
for u — 0.25, P — 2t^/v. Perturbations of a plane wave lead to defects only 
above the SN, whereas below the SN such perturbations have to be very large 
to overcome the saddle-type upper branch MAW. 

There are no SNs for parameters below the SN corresponding to P — > 00. Thus, 
starting from random initial conditions defects may only form at parameters 
above the SN of P — > 00. The SN of P — > 00 represents a lower bound for 
defect formation. 

For large systems the formation of defects depends on the local period of 
initial perturbations in a similar way as for v = [27,28]. The peak to peak 
distances of Lpx{x,t) are used to determine local periods p. In that context, 
the following explanation of defect formation has been proposed. Defects are 
observed in the phase turbulent regime whenever local structures, similar to 
MAWs, with spatial periods p larger then Psn occur in the system. Where Psn 
denotes the period of the MAWs at the SN (that coincides with the maximal 
MAW-period) for the considered choice of parameters. For larger values of u 
or C3 the SN occurs for smaller Psn as shown in Figs. 11 and 12. Therefore, at 
larger C3 local periods p beyond the SN and subsequent defect formation are 
more probable. In contrast to the case u — 0, there is only a short transient of 
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Fig. 13. (a) Bifurcation diagram as in Fig. 3 showing the minimum modulus of 
the MAWs. (b)-(e) Numerical simulations illustrate the dynamics near the SN cor- 
responding to the arrows in (a), (b) Plane wave perturbed at one point and (c) 
unstable saddle-type MAW plus noise converge to the stable MAW. (d) Unsta- 
ble saddle-type MAW plus a different realization of noise evolves to a defect that 
changes v to 0. (e) As (b) but beyond the SN which makes defect formation possible 
for arbitrarily small perturbations of the plane wave. Note the long living transient 
of a non-coherent modulation, (b)-(d) are at C3 = 0.5 below the SN and (e) belongs 
to C3 = 0.55 above the SN for u = 0.25, P = L = 2t:/v = 25.13. 

phase chaos in the simulations with nonzero initial Vi > vm- The distribution 

of local periods p of the perturbations is given by the realization of the noise 
in the initial condition. For local periods slightly above the relevant SN (as in 
Fig. 13e), the perturbation first increases to a modulation similar to MAWs 
and appears almost saturated for some transient time until finally a defect 
appears. This transient of defect formation becomes shorter as the distance to 
the SN grows. If initial conditions are prepared with z/j ^ z/^, then some local 
periods will be far beyond the corresponding SNs. The transients of defect 
formation are shorter in this case. 

These two observations suggest an interpretation of the curve u{t) represent- 
ing the temporal evolution of v during transients with z/j > Vm [20]. The larger 
the initial z/j is chosen the smaller is the final value vj . The time scales of 
competing processes have to be considered. Local defect formation will not 
instantaneously effect distant spatial locations along the system. Instead the 
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local change of the average phase gradient v via a defect will take a transient 
time to relax over the entire system. For z/j much larger than Vm-, defect for- 
mation happens on a short time scale and independently leads to defects at 
many different spatial locations before the relaxation of the decreased average 
phase gradient Vf <S vm can stop defect formation. For vi just above um de- 
fects form slowly and the reduced Vf can relax the phase gradient at distant 
locations before other defects occur. 

Let us now verify if the mechanisms proposed to explain defect formation in 
the phase chaos regime for solutions with z/ ^ still hold for the wound-up 
phase chaos regime. For z/ 0, defects form if and only if the period p of 
local structures is larger than Psn, where Psn is the period for which a SN 
occurs at the chosen coefficients Ci and C3 [27,28]. Wc have considered two 
cuts in the parameter space at C3 = 0.5 and C3 = 0.65 and investigated the 
distributions of periods p for solutions with average phase gradient v ~ vm- In 
particular, 50 realizations of a system of length L = 512 initialized as a plane 
wave with wavenumber v plus noise in the amplitude and in the phase have 
been considered. Then these different initial conditions have been followed for 
an integration time t = 500,000 — 1,000,000. The last part of the run has 
been examined in order to extract the length of the coherent structure with 
the maximal period Pmax occurring during the evolution. 

From these simulations we obtain the following: if defects occur then Pmax > 
Psn in all observed cases. However, it is not true that a defect is formed 
any time we observe Pmax > Psn- If we let the system relax for a long time 
(t = 500, 000) and measure Pmax, then the number of initial conditions leading 
to a Pmax > Psn without defect formation is noticeably reduced. For C3 = 0.65 
the maximal conserved phase gradient is um ~ 0.1 . In the simulations we do 
not observe a defect for u — 0.086 and u — 0.098 but in the first case only 2 % 
of all runs show Pmax > Psn, while in the latter case this percentage increases 
to 8 %. Increasing u the maximal period Pmax increases and more and more 
situations with Pmax > Psn are found upon approach to um- 

The differences to the z/ = case may be explained by the coexistence of 
chaotic and stable not-chaotic attractors. Depending on the initial conditions, 
the solution of the CGLE can evolve towards one or the other. Therefore 
the system may exhibit local structures similar to multi-hump MAWs that 
possess SN bifurcations at parameter values larger than those of single MAWs 
(compare Fig. 10). In that case, some periods p may even exceed Psn- 
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Fig. 14. Theoretical bounds for i^Mics) (filled circles) : SN of P ^ oo as lower bound 
(thick solid curve) and the presence of splitting at SN as upper bound (dotted curve) . 
Other curves as in Fig. 2. See text for details. 



3.2 Limit of wound-up phase chaos 



For random initial conditions with v in the narrow range between the SN of 
P — > oo and the existence hmit of MAWs (see Fig. 6), it depends on the 
specific realization of the noise whether a defect can form or a stable MAW 
results. No defects form below the line z/m(c3). In order to understand this 
observed limit z/m (cs) of wound-up phase chaos, it is sufficient to consider the 
SNs of single MAWs. 

Although initial conditions with large P beyond a SN could lead to defects, 
this is often prevented by the action of the splitting instability. Then the period 
is decreased before a defect can form. Following the SN curve in Figs. 11 and 
12 one encounters a transition from SNs with a splitting instability at large 
P to SNs without this instabihty at short P. Defect formation in wound-up 
phase chaos with v > 0.1 occurs for parameters where the splitting instability 
is not present near the SN, i.e. above the dotted line in Fig. 14. 

Fig. 14 summarizes the bounds found so far for the limit of wound-up phase 
chaos. The domain of stable plane waves at low C3 is limited by the Eckhaus 
instability (thin curve). Within the shaded area only supercritical HBs of 
different period P occur but no SNs. This area is limited by the lowest SN 
curve of P — > cxD (thick solid curve). No defects do form from random initial 
conditions within the shaded area. The dashed curve denotes the upper limit 
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of the existence domain of MAWs with any period. Saddle-node bifurcations 
exist in the window between this dashed curve and the thick sohd curve for the 
SN with P — > cxo. Defects will always form for any choice of initial conditions at 
parameters above the dashed curve. The dotted curve gives the transition from 
active (below) to inactive (above) splitting modes at the SN. This transition is 
computed by linear stability analysis along cuts like Figs. 11 and 12. Splitting 
inhibits defect formation below this dotted curve. Filled circles correspond 
to I'm obtained from numerical simulations as in [21]. Diamonds refer to the 
transition from chaotic (below) to non-chaotic (above) asymptotic states. 

The SN for P ^ oo (thick solid) is a lower bound for defect formation which 
also holds in the limit z/ = 0. The point Li marks the transition from phase 
to defect chaos studied earlier [27,28]. 

As long as the dynamics is regular (z/ > 0.1), the upper bound for the limit 
(onset of splitting at SN, dotted curve) of wound-up phase chaos reproduces 
well the numerical observations. For chaotic states with u < 0.1 defect for- 
mation eventually becomes possible even despite the presence of the splitting 
instability. This coincides with the increasing weight of the instability of single 
MAWs to attractive interaction of subsequent amphtude peaks. 



4 Experimental observations 

A variety of experimental observations in quasi one-dimensional geometries 
can be well interpreted by MAWs. These systems shall also serve for testing 
further properties of MAWs [37] . 

(i) For the oscillatory instability of Rayleigh-Benard convection patterns in 
an annular cell, B. Janiaud et al. report long living transients of modulated 
waves that eventually cause defects [4]. The underlying Eckhaus instability 
was found to be subcritical. Then, we expect the modulation to grow roughly 
exponentially whereas the long transient of an almost saturated modulation is 
similar to the dynamics near the saddle-node bifurcation as shown in Fig. 13e. 
In the latter case the Eckhaus instability is supercritical and stable MAWs 
may exist for nearby values of the experimental parameters, see also the dis- 
cussion in [38]. In all cases a single pulse-shaped modulation with the period 
equal to the cell length was present in the system and led to a defect. The 
excitation of several modulations per cell and thereby smaller period can sta- 
bilize the modulated pattern and provide more examples of the dynamics near 
the saddle-node bifurcation. 

{ii) Hydrothermal waves have been studied in ring-shaped cells [6] as well as 
in linear cells [5,7,8]. N. Mukolobwiez et al. report a supercritical Eckhaus 
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instability, stable MAWs with the period equal to the cell length and defect 
formation after a parameter change [6]. N. Garnier et al. observe modulated 
waves with both the wavenumber and the period of modulations being selected 
by one of the longitudinal boundaries [7] . Thermal or mechanical forcing at the 
boundary may yield more insight into the multistability of the two-parameter 
family of MAWs. 

[iii) In rotating Rayleigh-Benard convection Y. Liu et al. observe the subcrit- 
ical Eckhaus instabihty of a travehng wave sidewall mode [10]. The authors 
suggest higher order corrections to the CGLE in order to explain the observed 
discrepancy between the linear group velocity and the observed velocity of 
finite wavelength perturbations. However, this difference already follows from 
the linear analysis of the Eckhaus instability [16] and can be used as a test for 
the assumed coefficients of the CGLE. 

{iv) For the Taylor-Dean system 1. Mutabazi et al. report a stable MAW that 
they called "triplet state" because of the length scale ratio of modulation and 
underlying wave [11]. They also observe the formation of defects. Clearly the 
triplet state is just one realization of the two-parameter family of MAWs. 

{v) Finally we mention the heated wire convection [9], internal waves [12] and 
the oscillatory variant of the Belousov-Zhabotinsky reaction [13] where further 
investigations of the observed Eckhaus instabilities may also reveal MAWs. 



5 Discussion 

The bifurcation analysis of modulated amplitude waves (MAWs) has been ex- 
tended to nonzero average phase gradient {u ^ 0). Small amplitude MAWs 
("lower branch" -MAWs) of specific spatial period P exist between a super- 
critical Hopf bifurcation (HB) and a saddle-node (SN) bifurcation. The HB 
asymptotically reaches the Eckhaus instability from above as P goes to infin- 
ity. MAWs are a direct consequence of the Eckhaus instability of plane waves; 
they are obtained by a computer- assisted nonlinear analysis of this instability. 
We encounter SNs with decreasing values of P as Ci,C3 and v are increased. 
These SNs govern the formation of defects from random initial conditions 
as well as many aspects of the evolution of wound-up phase chaos. The SNs 
bound the existence region of MAWs in the Eckhaus unstable regime. 

A hnear stabihty analysis of MAWs revealed that they can be linearly stable 
even in systems of infinite size. These domains are limited by the interaction 
instability at low and the splitting instability at high values of the spatial pe- 
riod P of the MAW. The competition of the two instabilities drives wound-up 
phase chaos and determines the degree of chaoticity of the dynamics. For fixed 
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coefficients Ci and C3, the SN associated to P ^ 00 occurs at the lowest value 
of V. This SN establishes a "lower" bound for defect formation and thereby 
for the limit vm of wound-up phase chaos. The splitting instability can inhibit 
defect formation if the SN occurs at large P hence defects are created more 
frequently for parameters above a second curve where the splitting instability 
vanishes at the SN. Thereby an "upper" bound for the limit vm of wound-up 
phase chaos is obtained. Earlier numerical observations on vm{c^) are well 
reproduced for u > 0.1, respectively small C3. For u < 0.1, the description 
of phase chaos rehes on considerations similar to those already discussed in 
the hmit case u — [27,28]. Finally, several experimental observations were 
interpreted in terms of MAWs. 

We acknowledge fruitful discussions with M. van Hecke, N. Garnier, F. Davi- 
aud, W. van de Water and L. Kramer. AT was partially supported by the 
MURST-COFINOO grant "Chaos and locahzation in classical and quantum 
mechanics" . 
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